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ABSTRACT 


Long term numerical integrations were performed with a two-level 
model utilizing the quasi-geostrophic equation set. Friction was in- 
corporated through the vertical derivative of the eddy stress and heat- 
_ing was applied as a linear function of latitude. The long-term 
interactions between one or two waves in the zonal direction and the 
mean flow were examined. 

The stability of the initial mean flow was investigated with lin- 
earized equations and the most unstable wave numbers were determined. 
The nonlinear equations were integrated using the most unstable wave 
in the initial conditions. 300 day forecasts were made with various 
values of the heating, PB, wall separation, and wave number. A second 
wave number was introduced into the experiments in order to determine 
the most likely wave number for the fully evolved mean flow. Constant 
amplitude, baroclinic, propagating disturbances were obtained in every 
case where the wall separation was 4,000 km. The experiments with the 
wall separation of 8,000 km produced non-steady waves with fluctuating 


baroclinic and barotropic interactions. 
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I. INTRODUCTION 


In a classic experiment, Phillips (1956) performed a long-period 
numerical forecast with a two-level quasi-geostrophic model in a zonal- 
ly periodic channel on a Beta plane. Friction and heating were includ- 
ed in his model with the latter as a linear function of latitude. 
Although Phillips' experiment did not reach statistical equilibriun, 
he demonstrated that the principal features of the general circulation 
of the atmosphere could be reproduced with a simple numerical model. 
More refined studies involving the primitive equations have been car- 
ried out by Smagorinsky (1963), Mintz (1965), Manabe and Smagorinsky 
(1967), Manabe (1969), and Kasahara and Washington (1967) and (1971). 
These studies which more accurately model the actual atmosphere are 
sufficiently complicated so that it is difficult to understand the 
basic interactions between'the mean flow and the disturbances. 

For this reason, we wish to re-examine Phillips’ original simple 
model. In this study, Phillips' model is further simplified to repre- 
sent the x-variation of the disturbance with one or two Fourier re 
ents. This restriction is partially justified by the fact that one 
predominate wave number was observed. Similarily, in dish-pan experi- 
ments with inner cores, single wave numbers were observed in many in- 
stances; see Fultz (1959) and Hide (1953). Herring (1963) and (1964) 
has used the technique of examining the interaction between a disturb- 
ance of sinusoidal variation and a mean field in his study of convection 
between two rigid plates. His estimates of the heat fluxes between the 


plates were within 20% of the value obtained from laboratory experiments. 


This suggests that a simplified dynamic model of the atmosphere which 
involves the interaction ОЕ a simplified disturbance and a mean 
flow can give reasonable approximations of the northward fluxes of heat 
and momentum. 

Therefore, in this study, the complexities have been kept to a 
minimum to simplify interpretation of the interactions. The two-level 
model utilizes the quasi-geostrophic equation set. Fluid motion in the 
model is restricted to the mean flow and one or two wave numbers. The 
mean flow and Fourier amplitudes of the two disturbances are allowed to 
vary in time and y-space; hence an adequate description of barotropic 
as well as baroclinic interactions are possible. 

However, such a drastically simplified model must omit many details 
of the real atmosphere. In the case of this model, a number of physical 
properties of the earth's atmosphere are missing. Two waves in juxta- 
position will not interact directly, but can only influence each other 
through the modification of the mean flow. № oceans, moisture, moun- 
tains, or east-west heating gradient are allowed in the model. Conse- 
quently, this model appears to be one of the simplest models which can 


describe some of the principal mechanisms of the general circulation. 


= 


II. THE MODEL AND FORECAST EQUATIONS 


The model used in these experiments utilizes the simple two-level, 
quasi-geostrophic equation set. The atmosphere is divided into four 


layers of thickness Ap/2, numbered 0 to 4 in Figure 1. 
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Fig. 1. Two-level model. 


The earth's surface is assumed to be a smooth, flat, land surface 
with no oceans. Therefore, the vertical motion at the top and bottom 


of the atmosphere is taken to be zero. 


Internal friction is represented by the vertical derivative of the 


eddy stress T. 


The quasi-geostrophic vorticity and thermodynamic equations are 


9 : до _ : OT 
so + Y У + Вох - Е 5р = -М Ужв 51 
3 20 Й = |. AS 
Sr SA V DU Dc E 


where A - R/C, › S = net heating from radiation, and O, the static 


(2.1) 


(2.2) 


stability, is 


2= 2 
— “Ret RT dT 2 -2 -2 
с = m Var = E ео о-в cb (2.3) 
P8 Рв Р 


T is taken from a standard atmosphere. 


Apply Equation (2.1) at levels 1 and 3, which yields 


961 v | 5702 “Bik 7 | 
ot М ° 61 B Yi 5 > Др Е Ар! E ux Ше?) с 
Ot Be ©з Ex 3 E Ap za у ш” 5 12) 2. 


If we apply Equation (2.2) at level 2, we obtain 


besonders I 
Ba ae 50050 cue (2.6) 


The heating function, So» is restricted to a linear function in y 





as utilized by Phillips (1956). The model has heating in the south and 


cooling in the north. 


з, З НС > У (27) 


where w is the north-south wall separation distance. 

The amount of net heating is proportional to h which varies from 
0.001 to 0.016 kilojoules cont oan 

The eddy stress т, at the top of the atmosphere is taken to be zero. 


The eddy stress term T, at level 2, takes the form 


и . 
2 "88 ' р (2.8) 


where y is the eddy viscosity. If we follow Smagorinsky (1963), the 


finite difference form becomes 


т, - -(ру); (во,/4р) СМ - W) (2.9) 
where (0 ,8/4p) = 1/7.9 km! is the inverse depth of the 750-250mb layer 
50 | 
апа (07), ={225 (х 107% ton =“ a is an exchange coefficient. The 


500 
smallest estimates of (Pr), was proposed by Rossby and Montgomery (1935) 


for stable conditions, while the largest was utilized by Riehl (1951). 


The intermediate value, 225, was proposed by Palmen (1955). 


Thus, the stress term at level 2 becomes 
T, = -= с, NW - ND а (2.10) 


50 -6 -1 
where C, - 1225) x 10 /79 ton s (2-11) 
2 (soo 


The majority of the experiments were conducted using Palmen's exchange 
coefficient of 225. | 
The eddy stress term Т, ас the surface, level 4, as formulated by 
Phillips (1956) is assumed to be in the direction of the surface wind 
and proportional to the square of the wind speed. 
gs zb с Iw, I. (2.12) 
where Cy is the drag coefficient = 0.015. 


The surface geostrophic vorticity is approximated by the 750mb 


vorticity 


гор У ЧЕ 


о 
thus M = W 5 
So, by defining C, = 5 C Iv, | (2.13) 
PM (2.14) 


The average value of the 750mb wind speed was taken to be 5 m s 


Thus, our assumptions make Т, a linear function of wind speed with 


С, = 3.69 х 1073 ton s^! (2.15) 


Therefore, as boundary conditions at the top, internal, and bottom 


layers of the atmosphere 


о = 0, = 0 
т =0 

о 

T = 


25-6, 0%- М) 
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For convenience, the following quantities are defined 


+ 
= S033 
dy "ENS > | (2.16) 
+ 
® ЗЕТЕ G 
$. ж (20195 


implies 


= 


+ 6 %-%-% 


61 + 63 61 7 63 
ы онна ш а 
G1 = Sy + Sp G3 = Sy ~ Sr 


the mean quantities may be interpreted as applying at level 2. 
Integrate the hydrostatic equation d. - Q from level 1 to 3, 


which yields 
TERT 
6, "LM UP) | (2.18) 


Equation (2.18) shows that 0, is proportional to the mean temperature 


of the layer. 


< 


ТЕ (2.16) and (2.17) are introduced into the thermodynamic equation 


(2.6), we obtain 


w,OAp NS 
же НО ee о 


The mean vorticity equation is obtained by adding Equation (2.4) 





and (2.5) and dividing by two. 
9 ag 
Se Cur Ny Vout Wo E, + В Ум ^ 2% ко Ух т, (2.20) 
The thermal vorticity equation is obtained by subtracting Equation 


(2.4) from (2.5) and dividing by two. 
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9 е e — = 
Seen + N ge + We VO + ымы, - Е = 


ee (2421) 


2 24 4 


Solve (2.19) for о, and apply this solution to the thermal vorti- 


2 
city equation (2.21). 
Substitute the following geostrophic relations 


ya HIkxv $, се + у? ф 
о о 


into Equation (2.20) and multiply by f, ; which becomes the Mean Vorti- 


city Equation: 


vob, + = (ik x vào . ув i g Ok x 70-797) 


9%, 

2 2 = 
+В озу - kK, (7b, - 7) = о (2222) 
Introduce the geostrophic relations for V and C into (2.21) and 


multiply by f, » Which gives the Thermal Vorticity Equation: 


La ub e 1 akex vip ve? - й 
о 


ды 
e$ x vip * veo - 8, s * K (V^, - “фр 





AS 
зи _ 2 
tM 5 -REV б, (2.23) 
where 
um = 
bs = 9- = 3.48 x 1022 as (2.24) 
(Ap) 
С, 5 
= 4° 2 -6 -1 
к, 245 s 1.12 x 107" < (2.25) 
20,2 
=_2° = -6 -1 
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It can be seen that the frictional coefficient at the surface is 
roughly three times as large as the internal coefficient. 
The mean and thermal vorticity equations (2.22) and (2.23) are the 


prediction equations in the two unknowns 0, and b, . At this point, 


the geopotential height fields are restricted as follows: 


4, 


E(y,t) + A(y,t) cos kx + B(y,t) sin kx 


+ G(y,t) cos mx + H(y,t) sin mx (2227) 


br 


F(y,t) + C(y,t) cos kx + D(y,t) sin kx 

+ P(y,t) cos mx + Q(y,t) sin mx (2.28) 
where E(y,t) and F(y,t) are the amplitudes of the zonal mean flow; 
A(y,t), B(y,t), C(y,t), and D(y,t) are the Fourier amplitudes of the 
disturbance of wave number k, and G(y,t), H(y,t), P(y,t), and Q(y,t) 
are the Fourier amplitudes of the disturbance of wave number m. Both 
waves are in the 61 only. 

Since E and F are functions of y, barotropic instability is possible 


in this system as well as baroclinic instability. If the definitions 


of ф, апа 4, are substituted into М = tk х уф » we obtain: 
аб, o y 190. 
"M і, ду М і, Ox 
- 1 OE - 1 ӘЕ 
= ~ = — = = — = 2. 0 
Чм ЕЕ ду Ум i = (2224) 
581 94, ure 54, (2.31) 
T і, ду T і, Ox 
- _ 1 OF - | OF | 
un == Е ЗУ Ут = i T 0 (2.32) 
о 


where the bar C) indicates the x-average. 
At the north and south boundaries, the following conditions are 


imposed initially and for all time: 
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оао at y = 0, y=W 
өне ИЕ ——" 
т ee at y = 0, y= 


where W is the northern boundary. 

Substitute (2.27) and (2.28) into the mean vorticity (2.22) and 
thermal vorticity (2.23) equations. Now if we separate the various 
terms, and neglect all terms of wave number 2k and 2m and combinations 
of m and k, we obtain ten equations in ten unknowns. 


When we equate the coefficients of the cos kx terms of (2.22), we 


obtain: E 
2 2 2 3 3 
о 2), - E OBOE, ODOF _ ,O E SE Q2(59E , 9E 
Қт із шоу саз a 
En 2 
- В ВК - шо - К°) (А - С) (2.33) 
y . 
The sin kx terms of (2.22) give: 
2 2 2 3 3 
д 2 к . ДАДЕ O°C OF ОЕ, ОЕ 2 дЕ | (OF 
wer a Sy А3 A )] 
xX А ay2 ЗУ 57 ду ay! ay? e 
э 2 
+ B Ak - Әсер, - К/) (B ~ D) (2.34) 
The cos kx terms of (2.23) give: 
2 2 2 
ә 9 2 2 k ОЕ "р , OF 9°B OF 
= (—_-k“- a ту = S - == = 
TN, НОС Е УЗ “oy 52 5 к SEB (ke и?) 
3 2 
В 9 Fs] . B Dk + к, (9-02) (А-С) 
ду ду і ду 
2 
+К, (ck? - - (2.35) 
y 
The sin kx terms of (2.23) give 
2 2 
ә о 2742 k JE IC dE DA ОЕ OF 
Sega HD gE Sy O ан 
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The 


The 


The 


The 


E. ӘЗЕ o? 32p 


+ So + SEA] B. Ck + K, (—,-k) (B-D) + K, (Dk°- 22) (2.36) 
ду зу 2 E а 
terms independent of x in (2.22) give: 
2 2 2 2 
д 9 кво АВ QNA . 3 B о 
GE "zs [AR Баро | Со“ то 
dt ay) 2% dy - x^ wa) a? 
2 2 2 2 
cC T LI 
о ду ду ду ду 
OE OF 
E, S, Gs = (2.37) 


cos mx terms of (2.22) give: 


поо m (Ә2н ӘБ , 920 ӘЕ _ ЭЗЕ _ OF 2 


OE OF 
(—, - m )G = [— > Se Н--- Отто (HT + 9-21 


52 2 | 
- B Hm - К, С - що) (С - Р) (2.38) 
о ду 


sin mx terms of (2.22) give: 


2 2 2 3 
9 2 ш, Е ОР, ОЕ. ОЕ. 2 06) ЈЕ 
Зе m, ы. gn (Al 
о a y y 
+ В са - к,(9—- 2) (н - 9 (2.39) 


cos mx terms of (2.23) give: 
2 


3 32 2 2 т dE 379 | а ОБ 2. 2 ОЕ? 2 
(8—- п2- р2ур = m [SE - SQ u^- H(m - 4) 
dt з f. ду y ay? "ду 2 S y ду 
ЭЗЕ 2 32p 


3 F 3? 2 
233072 gil + PQ KE. ИОС) е КО Еше 5) (2.40) 
ду ду ду ду 


sin mx terms of (2.23) give: 


9 37 2.2. m, dp oF dfe , 3E ЈЕ 
wi = UO = Hu 7 ae 552 + SP nt зи Ae Sem” 18 9 
3 3 
+ E + Е E P Pm E K, п m )(H-Q)+ K, (qn^- 2^0 (2.41) 
ду ду ду ду 


15 


The coefficients independent of x in (2.23) give: 


<- 


2 2 2 2 2 
re a - m 
ду Өй” >. ду ду ду | 
2 
m à 9 e д Р de 36 2 
ze IE pS - 025 + u (HP - 69)] 
ош Su зу? ay ay бу 
AS, 2 
K, СЕ - 85 443.427 (2.42) 
y ду 


Equations (2.33) - (2.42) constitute the prediction formulas. 
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III. PROCEDURE AND GOALS 


The north-south grid distance Ay was set at 200km. A time step, 
At, of 1/2 hour was chosen to allow wind speeds up to 100 m goi with 
computational stability. 

Centered time differences were used for all quantities except those 
involving friction. The frictional terms were computed at time t - At. 
The first time step was a forward step. The Euler-backward time step 
was introduced every 24 hours to selectively damp high frequency waves. 
These waves are produced by the separation of solutions at even and odd 
time steps caused by the leap-frog time scheme. 

The finite difference scheme used is illustrated below with the 
function F, 


OF 1 


Sy | Зубні 7 8-1) 

2 

ӘСЕ 1 

— = — F. - 2F. Br 

552 (лу)? са о 

3 

OF 1 

S5 = — wo” 2F.,,.+F,)-(F, - 2, , +F, 
дуз 2 pa Fint pC} 2-1 ie 


where j is the grid index and Ay the distance between grid points 
(200km) . 

Equations (2.33) through (2.42) except (2.37) were solved by a 
Gauss pivotal elimination matrix method described by Richtmyer (1957). 
Equation (2.37) was solved by a direct marching process). 

The model was run as a Beta-plane experiment with periodic east- 
west motion between north-south walls of separation W, 4,000 or 8,000 


km. 


17 


The 8,000 km experiments were likened to the atmosphere. It was 
hoped the mean — MÀ from the north to south boundary 
at level 2 would be comparable to that of the atmosphere at 500 mb. 
The source of available potential energy, (S*P), should yield a value 
not much less than 60% of the atmosphere's total zonal generation of 
56 x 107 ^x3 52 7 formalized by Dutton and Johnson (1967). It should 
be-noted that this model contains no east-west heating gradient; there- 
fore, conversion of heat into available potential energy (APE) will be 
less than observed. 


Lastly, the winds at 250 mb should exhibit a jet structure and the 


750 mb winds portray the belts of zonal easterlies and westerlies. 
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ТУ. VERTICAL MOTION EQUATIONS 


Vertical motion at 500 mb is calculated through the thermodynamic 
equation (2.19). 


E 2 SC с. -0Е 2495 3 Ор пи OE A OE 
2, = Apo 2 Кх [5 - ғ Poy Еа + sin kx [37 + f Ke | 


+ сов шкід 7 є зви i12 * sin пи [99 4 Е оу Ма у?! 
о 
E 1 дА, дс OB др 
+ St” 2 | Дрої, | к у и 550 т = ду 
OG OP OH 90. | 
+ m[- м у У! | (4.1) 


all terms of wave number 2k or 2m and combinations of m and k have been 
neglected. 
The vertical motion, w, can be written as 
0=0+0 
where the average vertical motion is 


СОЕ a a 
® = Too Е Bet ae № 5,066 - 40) e» S ON - 09] (4.2) 


again, the barred quantity means x-averaged. 


The disturbance vertical motion, w', is given by 


w' =o! (cos kx) +o! (sin kx) +o! (cos mx) + w! (sin me) (4.3) 
where 

ш! (сов кх) = ES | cos [С - EOS : 2i (4.4) 
a (sin kx) = = {sin ta [$2 n ES В S51 | (4.5) 
w'(cos шк) = є | соз me (SE Е А + = $51] (4.6) 
0! (sin mx) = CS | sin mx [SQ + Tagi д En | (4.7) 


V. COMPUTATION OF AVERAGE NORTH-SOUTH WIND 


The average north-south flow computations coupled with the average 
vertical motion provides one with a fast method of calculating the cir- 
culation in the meridional plane; i.e. Hadley or Ferrel cells. 


If.we apply the continuity equation at level 1, we obtain 


до, _ 

Эр. = - у • м (5.1) 
dv D 

= ЕЕ Ls (5.2) 


Thus, the finite difference form is 


CI O 1 


а, ы 
Wp OD; 7% [ 5 ] (5.3) 


where (v); = 0) 
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VI. ENERGY AND ENERGY TRANSFORMATION EQUATIONS 


The energy equations were obtained from Phillips (1956). 


Mean Potential Energy per unit mass, 


2 
LZ. u й 2 


swf 
о 


Disturbance Potential Energy per unit mass, 


2 
= [&- & 
P = (ф b,) “dy 
вне 7 1 3 


Mean Kinetic Energy per unit mass, 


ч! 


El =2 , 2 
== — + 
4W je 75 
Disturbance Kinetic Energy per unit mass, 


A то 2 
К' === lib) + (v6) ] dy 
чт?) 1 3 


(6.1) 


(6.2) 


(6.3) 


(6.4) 


If we substitute the %, апа A, into the above energy equations, 


we obtain 








2 N 
Р —, > 
24 1-0 J 
2 N 2 N 
Р' (К) E, x (02 + 92) Р'(п) = ES 5 (Р? + 
ane“ 140 I 3 ane“ 3=0 J 
N 
- 1 ОЕ, 2 OF, 2 
K = A G9 
NES УО зв 
М 
EN! 27) 2 ода ов 2 902 
K' (k) m 2 pa Re (АЗ + Bit ar p^ > + ЗУ у. о у) б (— 
о 
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2 
0.) 


OD 27] 


М 
1 2 2 OG 2. ЭН. 2 
' EE жие > 
K' (m) "E T ІНЕ G + | + Р, + а; ) + (55 j te 
ӘР,2 _ 902 
E t Gp; 


where the subscripts k and m stand for wave number, and N is the 
number of grid increments such that 
= МДу 
Following the procedure of Phillips (1956), the energy flow diagram 


of Figure 2 should make the Ра Pow a of the energy trans- 


и 
оГ 


| 
(K’-€2) 


formation equations clear. 


Fee) 








Fig. 2. Energy Flow Diagram 
The flow of energy is in the direction of the arrows if the asso- 


ciated transformation ( * ) is positive. 


The energy transformation equations are: 
E. u?x N 
(S*P) = z © (S, °F.) 

ә 3903 J 





(Р.Р!) (К) = ше ` Аер = вес), * GE F 
әне 7 і- 1 


-0 

N 

(P'*K')(k) = ale х [C.w' (cos kx) + р fs no kx) J 
2/p N j=0 Jc] 
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М 2 2 2 2 
к) = ea ee ру 
25. М 10 yJ ду ду ду ду J 
2 2 2 2 
* G5, qo. 400 + А. В) | 
Oe CS ese 
— — -1 N == 
(РА) = + cr @,-a,) 
ApN ғ j 
-K N 2 
= 2 ӘСЕ 
(Кет, )= DM SUE 
2 М j=0 24 no] 
K N 2 2 
(K'e) (k) = —2 E Dc? 55, - (0—6 + 12-0) I 
2Nf j=0 J Oy ду J 
к -K N 2 2 
(ет, = En E, 
N і, j=0 J Oy Oy 
K N 2 2 
осот) ® ar ao), 
2МЕ, | 370 J эу dy” J 
2 2 
- s [ÊZ - 2D) -~ к(в-р)] 
J ду Эу Ј J 


The energy transformations of the m wave disturbance have been 
omitted for brevity. By replacing k, A, B, C, and D with m, G, H, P, 


and Q, the entire transformation set will be complete. 
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VII. INITIAL CONDITIONS 


A natural initial state would be one of no motion and constant tem- 
perature. Then the system could run until a steady state is reached or 
some type of instability occurs. However, to save computer time, the 
initial state was chosen to be near a critical state for development of 
the most unstable wave. 

When the mean flow model was allowed to run with no disturbance, 
the wind field reached a steady state; the Reynolds stresses balanced 
the heat input. This steady state was comprised of a meridional Hadley 
cell superimposed on the zonal flow. This Hadley cell mean flow was 
utilized as the initial condition. 

If we set > = 0 E A-D, с-0 - 0, 


д ДЕ ДЕ 
Ка dy y ay)” 


the F equation yields 0 


2 
the E equation yields 0 = ^m - y)/W - K, oF . 
: ду 
The solution for F is: 
уи уг у? 
каш K, ое 0 (7.1) 


This solution requires that E = F which implies that à, = ф, . In this 
case 


u. =0 andu. = 20 = Zu, ; 


Pedlosky (1964) derived a simple stability criterion for the mean 
flow. If the potential vorticity gradient 99) changes sign in either 


the y-direction or the vertical instability is possible. 
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For this model dq/dy is given by 








2— 5 
Su 2 
да! В 1 Ы ,- == - 
= - + —(u - о.) (722) 
9 о 2 2 3 
y ду 
sid ә, а-н 
ee - (ч, ш) (7.3) 
о 2 2 1 3 
У ду 
Although ue = 0 initially, easterlies and westerlies develop with 
time as observed in the earth's atmosphere. However, because Пра 5 
we shall neglect u, for this discussion. 
If we neglect и. in equation (7.3), we see that 
2 
= ge 
Ро 2E 


1 

is the critical value of the mean flow. Above that value, instability 
may occur for some wave number. If we solve for ч, ВИСТА value of 
A] may be derived through the Hadley cell initial condition (7.1). 


In this paper 


b ritical = 0.0011 


The disturbance initial conditions are 
A = 200 sin Ty/w 
в, C, D = 0 
G = 200 sin Ty/w after 100 days 


H, P, Q =0 
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VIII. RESULTS 


A. HIERARCHY OF UNSTABLE WAVE NUMBERS 

The objective of the first experiment was to find a hierarchy of 
wave numbers in order of wave growth for a particular mean flow. The 
mean flow was given by the Hadley cell solution with h = 0.004 kj tone 
E and it was held constant with time. Since the mean flow did not 
vary with time, the equations for the disturbance became linear. Each 
case was run until expotential growth was achieved. Figure 3 shows the 
graphs of wave number vs wave growth rate for two sets of conditions: 
B, - 0 and B, = 1.67 x ШЕТ 


growth rate shown is 
Amp (t,) 
Amp(t,) 


7 





growth rate = ln [ ] x 10 


| 


where Amp(t) = amplitude of the disturbance at time t. 
Pedlosky (1964) has computed the growth rate for this wind profile 
using different values for various parameters. 


Figure 3 exhibits how the Beta term stabilizes the longer waves. 


B. INTERACTIONS BETWEEN WAVE AND MEAN FLOW 

The second set of experiments were devised to test the growth of 
the most unstable amplifying wave, allowing interactions with the mean 
flow, as a function of heating. All experiments were conducted using 
Palmen's exchange coefficient of 225 x 1074 ton A а and 
heat = horitical (0.001 kj fous: qs as the initial condition for the 


Hadley cell. 


* 
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Each experiment was run for 300 days. The k wave was introduced 
immediately; the m wave (if existed) was introduced. at 100 days. 

The following tables were set up to briefly describe the character- 
istics of each wave number by way of symbolic nomenclature. The sym- 
bols are: 


D = wave dies immediately 


N = a.non-steady wave, disturbance energy 
oscillates up and down as the wave interacts 
with the mean flow. 

SW = a steady wave, total disturbance energy 
15 = ( ) d E 

Peak - the peak amount of disturbance 
energy = ( ) т^ 21% at ____day. 

E'(k) vs E'(m) - — — of k and m wave 
disturbance energy over the 200 day period 
when both waves are introduced. 

E'(m) > E'(k) - if che m wave disturbance energy 
did exceed that of the k wave, then the m wave 
disturbance energy, Е'(т) = ( ) fa ye 


on day. 


( ) kj n g^ is the average amount of conversion of heat 


с | 
e 
Fd 
! 


to APE usually computed after the mean flow has come to some 
degree of equilibrium with the disturbance. If the system has 
not run long enough to reach the equilibrium, then the value is 


given which best represents the average S*P thus far incurred. 
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C. u AND Up VERSUS HEATING FOR THE STEADY WAVE EXPERIMENTS 


The table below shows the dependence of the zonal flow upon the 


net heating rate for the Steady Wave experiments. 


heat (kj сой Т st ехрії uj sh u, (m 525) 
0.004 3 6.9 14.1 
0.008 4 7.1 15.1 
0.012 5 22 15.5 
0.016 6 452 16.0 


Table 5.11 


Table 5.11 lists the steady state values of the zonal flow. The 


rather small increases in Un with heat is best explained by the rapid 
4 


baroclinic response to the wind sheer decaying the unstable zonal flow. 
D. A BAROCLINIC, STEADY WAVE 


1.67 х 1071! 


4,000 km 
7 , exhibits a 
0.012 


о 


В 
м 
Experiment #5 | k 
h 


M 

А 

rt 
ии 


very interesting feature in this model: a steady, baroclinic wave. 
Figure 4 is a graph of mean kinetic energy, K, mean potential energy, 
P, mean energy, E, and total energy, E, for experiment #5. Initially, 
the mean flow grows rapidly, as shown on Figure 4, with negligible 


disturbance energy. At the 40th day, the zonal mean flow has reached 
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its peak and decreases rapidly while the disturbance begins to grow. 
By the 115th day, the лес atmosphere has reached a steady state. 
Initially, the maximum B = 4.9 m 87! апа чт = 2.4 ш se At the 
20th day, the maximum u, = 14.2 m 5375 (henceforth, whenever zonal mean 
flow is mentioned, it will mean the maximum speed in the north-south 


column) and m = 7.5 т 51, The 40th day exhibited the maximum mean 


flow energy, due to the maximum P, with u, = 21.7 m з апа u, - 11.1 


ms . However, the winds continued to increase until the 46th day, 
as attested to by the K graph, with En = 25.3 m un and e - 10.9 m d 

Figure 5 (E'(k)) portrays how the wave grew beginning at the 32nd 
day. By the 38th day, the wave responds to the strong baroclinic in- 
stability, as shown by the large up and grows expotentially.  Baro- 
tropic damping (K'*K) occurs exclusive of any barotropic growth the 
entire cycle of the wave, and exhibits large oscillations during the 
rapid growth of the wave. 

The wave continued to grow by baroclinic instability until the 56th 
day when the wave reached its peak energy with a disturbance energy, E', 
= 123 п s By this time, the thermal wind had decreased from a peak 
of 11.2 m 57 on the 42nd day to 7.5m PR on the 56th day, er 
the source of the baroclinic growth. 

On the 38th day, when the wave began its large baroclinic growth, 
barotropic damping begins to reach its greatest value as seen in Fig. >. 
Initially, as the barotropic damping oscillates about its mean of 


1.4x 10-5 


kj tons! the baroclinic instability is by far dominant. How- 
ever, by the 58th day, barotropic damping plus the loss of eddy kinetic 


energy, K', by frictional dissipation becomes larger than P'*K'; hence, 


baroclinic effects must be diminishing. 
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As the barotropic damping reaches steady state, it appears to come 
into equilibrium with бе. ай growth, hence, the steady state of the 
disturbance energy. 

Figure 6 is a graph of the conversion of net heating to available 
potential energy (S*P). The graph shows the rapid rise of S*P to a 


ы” ne st in a vertical column (1.6 x 107 "kj e 


peak of 16 x 10° 
ar days, »then its decrease to the steady state value of 8.8 
x 107 ^kj penal by the 104th day. 

At the end of the 300th day run, the entire atmosphere has been at 
a steady state equilibrium since the 115th day. The energy conversion 
terms feeding P, P', K, and K' were constant with respect to time, as 
were the energy terms; therefore, the conversion terms should balance 
for each energy term. Figure 7 is a pictorial display of the 300th day 
energy flow after Phillips (1956), the only conversion not balancing 
being the (кет,). 

In summary, the а. wave is achieved by the conversion of P to 


K' which is ultimately balanced by frictional losses and barotropic 


damping through the action of the Reynolds stresses. 


E. A NON-STEADY WAVE WITH BAROTROPIC EFFECTS 


В = 1.67 х 1071! 
w = 8,000 km 
Experiment #12 К = 6 , exhibits another 
heat = 0.016 


quite interesting feature, a wave initially amplified by baroclinicity 
which then grew by barotropic effects. 

Figure 8 is a graph of the mean flow energies and total energy for 
experiment #12. Notice the large and small oscillations in each energy 


graph and how the total energy graph is smoother than the mean flow 
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energy graph. This is due to the disturbance drawing energy from the 
mean flow and vice versa. © when the disturbance decreases, the mean 
flow increases and vice versa which tends to smooth out the total 
energy graph. 

On the 20th day, the wave grew rapidly by baroclinic effects. At 
that time, u, = 9.3 а с. uj = 18.1 п le das. 

m = 10,9 п si ug is = 30.0 м T Since only barotropic damping was 
occurring and the mean flow increased since the initial 10 days of 
growth, baroclinic effects amplified the wave to its initial peak of 
Eun = савеза days. 

Figure 9 is a graph of the disturbance energy and the barotropic 
instability. It may be seen that the barotropic growth term feeds or 
damps the disturbance. At day 30, the wave has peaked due to baro- 
clinic growth. From the 3lst to the 33rd day, the baroclinic instability 
subsided and the barotropic damping destroyed over half of the wave 
energy and increased the mean flow to ча = 15.3 м з and DS 38.2 m 
E 

The 34th is the first day that barotropic instability occurred and 
although it only lasted for two days, the disturbance peaked at 77 a 
E, On the 37th day, barotropic damping again diminished the disturb- 
ance until the 40th day. On the 37th day, E = 14.1 m and hy = 
37.9 m s while on the 40th day, ч = 17.2 т зві and m 42.6 m Su. 

To further supplement this discussion, the K graph is given on the 
bottom of Figure 8, without its energy scale. The graph helps to por- 


tray the barotropic damping effect, gain of mean flow energy to the 


disturbance and vice versa. 
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On the 40th day, the barotropic instability causes explosive ampli- 
fication of the wave. — in the 42nd and 43rd days, barotropic 
damping occurs and explosive baroclinic growth of approximately five 
times the average P'*K' takes over. 

This baroclinic growth can be seen on the small growth of the K 
graph from the initial loss due to barotropic growth from the 40th to 
“the 42nd day to the-small-gain from the 42nd to the 44th day. 

From the 44th to the 53rd day, growth is a combination of baro- 
tropic and baroclinic effects. Possibly the barotropic term is domin- 
ant, as with each barotropic damping, the disturbance graph changes 
slope or declines, hence a greater dependence on barotropic growth. 

The mean flow has declined greatly due to both a baroclinic and baro- 
tropic drain on it, as shown by the 44th day Up = 14.2 m se and 
5, = 40.5 m 571 while on the 53rd day, u,= 11.9 ms “and 3,5 23.9 m N 

From the 53rd to the 82nd day, the wave damped for the most part 
due to barotropic ее. although a small amount of baroclinic damp- 
ing does occur. 

The disturbance peaks on the 56th day, due to barotropic and baro- 
clinic growth, with the K graph reaching its minimum on the 55th day. 
Again, baroclinic effects do create growth in the wave as seen when 
barotropic damping is occurring and the disturbance pauses from its 
decay to grow for a day, then continue its decay. 

Figure 10 is included to observe the variation of S*P term in a 
non-steady wave. Figure 11 is included to show the many oscillations 
of the P*K term (thin line). An interesting phenomenon noted on Figure 


11 is the comparison of the P*K conversion with the K term (heavy line, 


without energy scale). The two terms tend to counteract each other 
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rather than reinforce each other. In the same fashion, K'*K also 
counteracts the P*K — with  K'*K on the average dominating 
P*K. Hence with barotropic damping, (кк > 0) the P*K conversion 
usually becomes negative. The reverse also occurs. 


Figure 12 is a graph of the u, flow with respect to time. The 


1 
graph begins when the disturbance and mean flow reach a state of 
-equilibrium after.the.initial .growth period of .the disturbance on the 
mean flow. The graph runs from 80 to 110 days. The easterlies present 


at 250 mb are attributed to the fact that no horizontal diffusion is 


present in the vorticity equation. 


F. A TWO WAVE EXPERIMENT 


= 1. х 1 
Ж 1.67 0 
Experiment #12 ү Е ы N ‚ 15 one of the 
m =7 
heat = 0.012 


few experiments conducted where E'(m) > E'(k). The explanation for the 
dominance of the m wave over the k wave, in opposition to the results 
of section V. A., is attributable to the en of the mean zonal 
flow by the k wave. Figure 13 is a graph of E'(k) and E'(m) versus 
time. Figure 14 is a graph of the K'*K conversions for each wave. It 
should be noted that, on the average, both conversions act in the same 
direction in quite varying degrees. 

The Pedlosky criteria (7.2 and 7.3) proved to be a necessary condi- 
tion for baroclinic and barotropic instability. The baroclinic criteria 
was met for the initial condition and throughout the forecast. The 
barotropic criteria was met continuously after the 34th day. At no time 


did barotropic growth occur when the Pedlosky criteria was not met. 
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Between the 162nd and 164th day, the E'(m) becomes larger than the 
E'(k) and remains such throughout the run. At this time, both waves 
are experiencing barotropic instability. However, the m wave is grow- 
ing and the k wave is decaying. 

The time frame, from 162 to 164 days, is explored further on 
Figures 15, 16 and 17. Figure 15 and 16 is a graph of P'*K' (bold line) 
and K'*K(thin line) for waves k and m respectively. A unique differ- 
ence is apparent between the two waves. The baroclinic conversion 
P'*K' never approaches zero with the m wave but does many times with 
the k wave, and even becomes negative. When the baroclinic term be- 
comes negative, baroclinic damping draws from K'. Usually, the P'*K' 
conversion was positive, thus the waves were experiencing baroclinic 
growth. On day 164, both waves are experiencing barotropic growth; 
however, the k wave is undergoing baroclinic decay and the m wave is 
experiencing baroclinic growth. The energy flow chart is shown below 
to help explain the growth of the m wave and decay of the k wave. 

| х | | [x ] 
| | 


1.41 Jr 
[ P J-ois-[ &' ]-038 Pa ]— 0.44 
i 
tn 052 
k wave m wave 
net change to K'(k) net change to K'(m) 
= +.61kj Poe su = +.92kj Sone = 


E  ——M + 


Possibly this difference in net inflow to the K' term for each wave 


gives the growth advantage to the m wave for a ten day period. 
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If we look at the waves on the 164th day, we find an interesting 
feature of the barotropic growth. Although the k wave K'*K is largest 
(negative), the m wave is experiencing the growth through baroclinic 
instability. Since the conversion terms are the meridional average, 
it is possible the barotropic growth of the m wave in a sector is sig- 
nificantly greater than barotropic growth of the k wave in a similar 
region. This possibility is explored in Figure 17, a graph of the à, 


phase angle of the k and m waves and u The wave is experiencing 


1 
barotropic growth if the tilt of the wave is opposite the shear of the 
mean flow. Figure 17 depicts the k wave growing barotropically south 
of the zonal maximum of u, to the easterly jet, and the m wave growing 
barotropically north of the zonal maximum to about y = 6,000 km, thus 


demonstrating the greater barotropic instability of the k wave. 
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Fig. 17. Phase Angle and v, vs y for Exp #26. 
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IX. CONCLUSIONS 


In order to facilitate physical interpretation of the results, the 
model was restricted to a mean flow and one or two waves. A steady 
state solution which was independent of x was obtained. This flow was 
"of the:Hadley type with zero mean velocity at level 3 and a parabolic 
wind profile at level 1. The stability of this mean flow was investi- 
gated by performing initial value integrations with the linearized 
frictionless equations. The calculations were continued until expon- 
ential wave growth was achieved. At this point, the growth rates for 
a range of wave numbers were calculated. The growth rate of the waves 
was determined as a function of wall separation and B. The growth rate 
was reduced by the B effect and increased greatly by large wall separa- 
tions. 

The nonlinear frictional equations were integrated from an initial 
State which consisted of the Hadley mean flow upon which was super- 
imposed the wave number of maximum growth rate. For each combination of 
heating, B, and wall separation, a 300 day prediction was made with the 
optimum wave number from the linear study previously mentioned. Each 
of these forecasts was repeated but at 100 days a new wave was intro- 
duced. The purpose here was to determine the largest amplitude wave 
which would exist on the mean flow, the latter having been modified by 
nonlinear interactions by the first wave. A total of 45 experiments 
were made as a function of wall separation, heating, В, and wave number. 

Constant amplitude propagating disturbances were obtained in every 


case where the wall separation was 4,000 km. These experiments 
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contained cases where B = 0 апа В # 0 as well as a large range in heat- 
ing rates. The steady waves gained energy through baroclinic conver- 
sion processes and lost energy through barotropic decay and friction. 

The remaining experiments with the wall separation of 8,000 km pro- 
duced non-steady waves with fluctuating baroclinic and barotropic in- 
teractions. The roles of the different types of instability in the 
growth and decay of the waves were. easily distinguished. When one in- 
stability was causing strong growth, the other instability was usually 
a minimum growth conversion or was decaying the disturbance. The baro- 
tropic and baroclinic conversions tended to be out of phase with one 
another, one a maximum - the other a minimum. The Р.К conversion, 
however, tended to be a constraint on the K term and adjusted to the 
barotropic term. With strong barotropic damping, the Р.К was negative, 
while with strong barotropic growth, the P*K conversion was positive. 
Hence with strong barotropic decay, the Ferrel cell became dominant over 
the Hadley cell. 

Over the long term mean, the waves obtained a net gain of energy 
through P'*K' and a net loss of energy through K'*K, The mean tempera- 
ture gradient from the south wall to the north wall and the S*P conver- 
sion were directly related to the heating and wall separation. With 
the largest heating and wall separation, the temperature gradients were 
somewhat greater than the north-south 500 mb temperature gradients in 
the earth's atmosphere. This is related to the larger than normal 
heating used to offset the lack of east-west heating gradient. Even 
with the largest heating and wall separation, the conversion of heat to 
available potential energy, S*P, is only about 65% of the average for 


our atmosphere of 56 x 10° "kj ШО QS as estimated by Dutton and 
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Johnson (1967). This underestimate is also the result of the lack of 
east-west heating. 

The zonal mean flow at level 1 exhibited a definite westerly jet 
in the mid-latitudes, which sometimes split into a double jet. Due to 
the lack of a horizontal diffusion term, the horizontal gradients were 
sometimes quite large and occasionally an easterly jet would appear for 
a short time. The zonal mean flow at level 3 exhibited the belts of 
zonal easterlies and westerlies typical of the earth's atmosphere. 

In this study, we found with a wall а Ва of 4,000 km only 
steady baroclinic waves, and with a wall separation of 8,000 km 
fluctuating baroclinic and barotropic interactions yield non-steady 
waves. It was determined that the predominate wave number in the fully 
developed flow was in every case close to the value determined by the 
linear stability analysis. In the actual atmosphere, where walls are 
not present, barotropic instability may be important in giving day to 


day fluctuations in the mean flow. 
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